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Abstract 

We show that a breakdown of the universahty of the gravitational cou- 
phngs to different neutrino flavors can be tested in long-baseline neutrino- 
oscillation experiments. In particular we have analyzed in detail a proposed 
experiment at SOUDAN 2 with beams from the Fermilab Main Injector. It 
turns out that we can study both masses of neutrinos and such a breakdown 
with sensitivity to the order of 10~^^ by investigating the energy spectrum 
of the resulting muons. 
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It has been pointed out^'^ that a breakdown of the universahty of the 
gravitational couphngs to different neutrino flavors could lead to neutrino 
oscillations. In particular the authors of Ref. 2 studied the possibility in 
which the solar neutrinos (see, e.g., Ref. 3) can be used to test this kind of 
breakdown of the universality. Since the flux of the solar neutrinos is rela- 
tively small and the energy spectrum is beyond our control, the utility of the 
solar neutrino for this purpose is limited. In this paper we propose a possible 
long-baseline experiments of neutrino oscillations to test the breakdown of 
the universality of the gravitational couplings to neutrinos. In case of the 
long-baseline experiments with an accelerator, we have a larger flux than that 
of solar neutrinos and in principle we can change the energy spectrum of the 
neutrino beams, and therefore more information, if any, can be obtained on 
neutrino oscillations. As we will see, the breakdown of the universality of the 
gravitational couplings to neutrinos of different flavors leads to a violation of 
Einstein's equivalence principle (see, e.g., Ref. 4) which states that all the 
laws of physics must take on their familiar special-relativistic forms in any 
and every local Lorentz frame, anywhere and any time in the universe. In our 
case, we show that we can probe the magnitude of the breakdown of Einstein's 
equivalence principle to the order of 10~^^, assuming that there are neutrino 
mixings. Among various experiments to test the equivalence principle (see 
e.g., Ref. 5 for a review), there have been few tests of Einstein's equivalence 
principle for neutrinos^. The universality of the gravitational couplings that 
we study in this paper is of different type from these experiments in the past, 
so our discussions here are complementary to them. 

In this paper we assume that there are two neutrino flavors which have 
different couplings to gravity and that the eigenstates of these different grav- 
itational couplings do not coincide with those of the electroweak flavors. 
Throughout the present discussions we consider neutrino oscillations between 
two flavors for simplicity. Let us start with the following Lagrangian of two 
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kinds of neutrinos 

+ e(G2)l72K^(G2)7a/^^(G2)-m2]z/2 

+ (interaction terms with electroweak gauge bosons), (1) 

where we have included mass terms to keep generality, e"^(G'j) {i = 1,2) 
are the vierbein fields of some background metric with different Newton con- 
stants Gi {i = 1,2), and e{Gi) = dete^(Gj). For simphcity we assume 
that the eigenstates of the gravitational couplings coincide with those of the 
masses. Notice that even if these neutrinos are massless, we cannot rotate 
these two fields so that these are the eigenstates of the electroweak theory, 
since the gravitational coupling terms are not invariant under the rotation 
in the fiavor space. Note also that each term in (|^) is consistent with lo- 
cal Lorentz invariance, general covariance, and the CPT invariance, as the 
gravitational couplings to particles and anti-particles are the same. Since 
the gravitational coupings for these two kinds of neutrinos are different, even 
if we choose a coordinate system in which the Dirac equation for ui in (|^) 
becomes the one in a flat space-time, the Dirac equation for z/2 in the same 
coordinate system does not necessarily do so. Thus Einstein's equivalence 
principle is violated in (|l]). The situation here is similar to that of the gauge 
theory where the gauge invariance is explicitly broken, and physics does de- 
pend on which gauge we choose. So we are forced to choose one particular 
coordinate system from which we start. The most natural choice in our case 
seems to be the coordinate system which is at rest on the Earth. This is 
because neutrinos which we observe are created from the accelerator and an- 
nihilated near the detector, and both equipments are flxed on the Earth. So 
we take the coordinate system which are moving together with the Earth, 
and we choose as our background the so-called interior Schwarzschild metric 
(see, e.g., chapter 11.7 in Ref. 7) whose curvature is entirely caused by the 
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gravitational field due to the Earth. We will discuss the issue of the choice 
of the coordinate systems again later. 

The configuration of the long-baseline experiment we will discuss is de- 
picted in Fig. 1, and the neutrino beams go underneath the ground along the 
geodesies. First let us consider the Dirac equation of left-handed neutrinos 
without any flavor in the interior Schwarzschild background: 



{ie''>'jaD^-m)^ = 0, (2) 
where e^^ is the vierbein of the interior Schwarzschild metric 

ds^ = {e^fde - {el)V - {elfde'' - (ej,)^ V (3) 
and is given by 



H3 

D^Tp = ((9^ — Itu^afcO""^)"^ is the covariant derivative acting on a spinor 
uj^ab is the spin connection given by ej'j.cij^j^ = 9[^e"], and a in (|) is the 
Schwarzschild radius. 

One typical dimensionless parameter in our case is ER, where E is the 
energy of the neutrino, and R is the radius of the Earth. For i?=10 GeV and 
i?=6,400 Km, ER ~ 3 x 10^^, and derivative terms in the spin connections 
are all of the order of 1/ ER, so we will neglect them throughout this paper. 
In this approximation the Dirac equation becomes 



= (5) 



Since we consider neutrinos in the ultra relativistic limit E ^ m, we have 
only to discuss (j^) along the geodesies for massless fields. The geodesies of 
the interior Schwarzschild metric is given by 

R^ COS^ 6.2 Set / 2 e- • 2 X 

= sm ^ + 77^ (cos d - sm ip) (6 j 
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to the order of a/R, where 5 is half of the angle ZAOB in Fig. 1. We 
solve d^) on the plane 9 = 7r/2, and we remove the time dependence by 
ipixjt) = e~*-^*x(2^)- Oil the geodesies 



oO\-l 



m + 



becomes 
^ r dip 



X 



0, 



(7) 



where w = dlnr{ip)/dip, r = r{ip) is defined in @, and we have taken only 
the tangential component into consideration. Following the convention of the 
Dirac matrices by Bjorken and Drell^, and denoting = (Xi? X2, Xs? Xi)y it 
is easy to show that (|^) can be rewritten as 



f Xi + iX2 
ir d(p\X3- 



1 + 



\ 1 + iei: 



-2^2 



Xi + iX2 
X3 - iXi, 



where U is a certain 2x2 matrix. It turns out that a derivative term 
r^^dU^^ / dip is of order 1/ER which is extremely small, and hence we get 



1 d 

ir dp 



V 
V 



-2^2 



1 + ^2 



-2y^2 



0"3 



V 
V 



(9) 



where (z/, v) = (xi + ^X2? Xz^^X-diU^^)^ . v and v correspond to the forward- 
going and the backward-going energy solutions, and we will consider only 
the forward-going solution in the following, a in eq. (|^) is the Schwarzschild 
radius of the Earth which is about 9 mm, so we expand to the first order 
in ajr. Thus, in the ultra relativistic limit E we obtain 

(10) 



Idv 
i dx 



B 



m a 



cos 6 — 2x cos 6 



where we have defined a variable x = tan(y9 — vr/2), and B = ER cos 6 is a 
very large number. Now let us go back to the Lagrangian (|lD with two kinds 
of neutrinos. This time we assume that the vierbein fields in the Lagrangian 
(|I]) are those of the interior Schwarzschild background (^. From the previous 
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discussion, it is straightforward to see that the Dirac equation for (|^) is given 
by 

1 



1 d 

i dx 



+ + (A + A)* 



4^2 



COS 5 — 2a; cos 5 



A(x)a3 
11) 



where 



A(x) = ^ + A/$ - cos^ 5 - 2a;2 cos^ 5] . 



(12) 



k A f * I 

Here Am^ = — is the difference of the masses, we have defined the 
Newton potential $ = —GM/R on the surface of the Earth, and we have 
also defined the difference A/ = /2 — /i of the dimensionless gravitational 
couplings of the two neutrino species 

ai/2R \ 
a2/2R ) 




GiM/R \ 
G2M/R I • 



(13) 



The equation ( |TTD can be easily integrated from x = — tan 5 to x = tan 5. 
Now let us introduce the flavor eigenstates Va-, i^b of the weak interaction 

by 

(14) 

Then the probability of detecting a different flavor Vb at a distance L after 
producing one neutrino flavor Va is given by 

L2 





P{va = sin^ 26 sin" 



f Am^ ^ A/$ 



EL 



(15) 



_V4E2 ' 2 y ' 6i?2 

where L = 2i?sin 6 is the distance AB in Fig. 1. This formula applies to the 
transition between z/^ and Ur, where no MSW effect^ is expected to occur. 

In case of the transition between and u^, we have to take the MSW 
effect^ into consideration, and the Dirac equation is modified as 



1 d 

i dx 



B 



A(x) cos2e 



A(x) sin 29 



V2E 



A(a;)sin 20 
-A(x) cos 2^ + ^ 



,(16) 
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where Gi? is the Fermi couphng constant, and N^. is the density of electrons 
in the Earth. (|16|) can be solved in the same way as before by introducing 
the variables 

GpN 

Aj^{x) cos29n- = A(x)cos2^^ -^=-^ 

A^(x)sin2^^ = A(x)sin2e. (17) 

Note that Qj^ does depend on the variable x in this case. It is easy to integrate 
(O), and we have the transition probability of detecting z/g at a distance L 



from the source of z/^ beams 



(/■tan (5 \ 
/ dx An{x) , (18) 
v/ — tan (5 / 

where we have used the fact Oj^{x = tan5) = 9j\f{x = — tanS) and Ajv, 
are defined through (|l^) . We have performed the integration in the exponent 
in (0) numerically. 

Here we would like to comment again on the dependence of our results on 
the choice of the coordinate systems. As we mentioned earher, equations (|^) 
and (0) depend on how we choose a coordinate system, since the Newton po- 
tential terms in (y) and (p!T| ) are changed if we switch to a coordinate system 
which moves with acceleration relative to the Earth. In Ref. 10 it was ar- 
gued that one could derive stronger bound on a breakdown of the equivalence 
principle by using the contribution to the Newton potential from the super- 
galactic cluster. However, this argument holds only when one assumes the 
Lagrangian (^ in the coordinate system which is at rest in the supergalactic 
cluster. Such a coordinate system is different from our coordinate system, 
i.e., a different choice of gauge is taken in Ref. 10. Throughout this paper we 
take the background (§), and we make our analyses below using this choice 
of the background. Since any other contribution to the Newton potential 



is negative definite, the absolute value of $ in ([15D with our ansatz is the 
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smallest among all possibilities. Hence our choice gives the most conservative 
bound on A/. 

Let us now consider a proposed long-baseline neutrino-oscillation ex- 
periment which will be performed with z/^ beams from the Fermilab Main 
Injector^^. Our discussions here are analogous to those by Bernstein and 
Park^^. Although there may be several factors which cause the systematic 
errors as has been emphasized in Ref. 12, we will not discuss this issue in this 
paper. We assume that the Fermilab Main Injector neutrino beams have an 
energy spectrum which is given by Fig. 6.27 in Ref. 11, but we extrapolate 
the graph in Ref. 11 naively for neutrinos of the energy larger than 70 GeV. 
The average of the energy of the neutrino beams is typically from 10 GeV to 
20 GeV. If Einstein's equivalence principle is violated, the higher the energy 
of the neutrino beams becomes, the more probability of neutrino oscillations 
we have (see eq. (p!5D), unlike in case of neutrino-oscillations due to masses. 
We choose a value of the distance L=800 Km which is motivated by the 
SOUDAN 2 experiment. We mainly consider the so-called disappearance 
experiments^^'^^, in which the initial z/^ flux at some short distance and the 

flux at the detector are measured by detecting muons which are created 
from charged-current interactions. The probability of detecting muons at a 
distance L is given by 1 — e(z/^ Ua) (a = e or r) where 

eivn Ua) = -^^^ — - — TTTT^ (a = 6 or r).(19) 

Here E and q are the energy of the incident and the outgoing muon, 
respectively, F{E) is the flux of neutrino beams, nxiq) is the effective number 
of target nucleons, we have modeled the detection efficiency function e(g) 
for muons with a step function, as in Ref. 12, and we have used the y 
distribution for deep-inelastic scattering of neutrinos to determine the cross 
section da{E,q)/dq of charged-current interaction^'^. P(z/^ Ve or ^t) in 
(p!9D is either P(z/^ — > z/e) or Pivp, Vr)-, depending on whether we consider 
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oscillations z/^ ^ z/g or z/^ ^ z/,-. 

We have studied the quantity e for various cases. In case of <-> Ur 
oscillations, there is no MSW effect, and we have obtained the results for 
e(z/^ Ur) from ([T5|) and (pj]). We have plotted, assuming no signal for beam 
attenuation up to e, the excluded region of the (sin^ 29, A//10~^^) plane for 
z/^ Ur oscillations with the detection threshold energy Eth=lO GeV of 
muons in Fig. 2. We have also calculated the excluded region for -Eth=5 and 
20 GeV, and the results are almost similar to the case of i?th=10 GeV. In 
Fig. 3 we have plotted the excluded region of the (sin^ 29, plane 
for z/^ ^ Ue with Eth=10 GeV. To get this region we have used ([T8|) and (|T9|) , 
where the MSW effect is taken into account. We have used a value for the 
density of electrons: = 8.05 x 10^^ electrons/cm^ which is constant on the 
entire trajectory of neutrino beams in the present case (see, e.g., Ref. 14). 
Comparing Fig. 2 with Fig. 3, we observe that there is a slight difference in 
'^^i ^ oscillations because of the MSW effect. We have shown the result 
with negative A/ in Fig. 3, because (p!7| ) shows that negative A/ gives larger 
|AAr(a;)| and therefore e(z/^ Ue) with negative A/ is larger than that with 
the same Am^ and positive A/. Notice that e(z/^ z/^) with Am^ = and 
—A/ is the same as e(z/^ z/^) with Am^ = and A/, as is obvious from 

One important feature is the dependence of e on the detection threshold 
energy E'th of muons. As we mentioned earlier, it is more advantageous 
to look at neutrinos (and therefore outgoing muons) of higher energy to 
investigate a breakdown A/ of the universality of the gravitational couplings. 
We have shown e(z/^ — > z/g) as a function of Eth in Figs. 4 and 5 with sin^ 29 = 
0.5 and with Am^ = 1 x 10~^eV^,0.5 x 10~^eV~^, and 0. The parameters 
are A/ = for Fig. 4 and A/ = —0.5 x 10~^^ for Fig. 5, respectively, 
^(t'/i J^e) with A/ = in Fig. 4 decreases as Eth increases, since the 
argument of the phase (=the mass squared) in case of ordinary MSW effect 
(see (0)) is suppressed by the neutrino energy E^, (oc l/E^). Fig. 5 shows. 
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on the other hand, that e(z/^ Ve) in the presence of A/(= —0.5 x 10"^^) has 
a conspicuous difference from the case with A/ = 0. Thus, if we look at the 
i?th-dependence of e(z/^ Ue), we can determine Af with higher accuracy. 
Because of the interaction term GpN^. with matter, e(z/^ Vg) is in general 
larger than e(z/^ Ut), but the qualitative features for the results of ^ z/^ 
oscillations are the same as those of ^ v^- The dependence of e on i?th is 
also very useful, if we can change the energy of the neutrino beams. Namely, 
by comparing the results with different two values of energy of the neutrino 
beams, we could establish the existence of non-zero value of A/. If the 
statistics is good enough, then we could even measure masses of neutrinos, 
by subtracting the effects of A/ ^ from the data. This is one advantage 
that long-baseline experiments have over solar neutrino experiments. 

Another important point is the expected numbers of observed muons. In 
case of the planned experiments at SOUDAN 2, we estimate the numbers to 
be approximately 16,000 events/year for £'th = 10 GeV and 6,000 events/year 
for E^\^ = 20 GeV, which are significantly larger than those of solar neutrino 
experiments where typical numbers are several hundreds events/year. So also 
in this aspect long-baseline experiments are promising. 

In this paper we have proposed long-baseline experiments to test the 
universality of the gravitational couplings of neutrinos, and we found that 
we could probe the dimensionless parameter A/ as small as 10~^^ which is 
smaller by a few orders of magnitudes than the upper limit on a breakdown 
of the equivalence principle from different types of experiments. Although 
we have not evaluated systematic errors in detail, we hope our analysis will 
stimulate and motivate long-baseline experiments in the near future. 

Noted Added 

Toward the completion of our paper, we became aware of the work by 
Pantaleone, Halprin and Leung^^, where the similar topics has been discussed 
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from a slightly different viewpoint. 
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Figures 

1. The cross section of the Earth. The accelerator is located at point A, 
and the SOUDAN 2 detector is located at point B. The distance L 
between the points A and B is about 800 Km, and the radius R of 
the Earth is about 6,400 Km. The curve from A to B in the Earth is 
geodesies in the background of the interior Schwarzschild metric. 

2. The excluded region in the (sin^ 2^, A//10~^^) plane for <-> Uj. 
with e=l% (solid), 3% (dashed) and 10% (dotted hne). The detec- 
tion threshold energy of muons is 10 GeV. The upper and right side of 
the curves is excluded. 

3. The excluded region in the {sin^ 26, —Af /lO^^"^) plane for z/^ ^ 
with e=l% (solid), 3% (dashed) and 10% (dotted line). The detection 
threshold energy of muons is 10 GeV. The upper and right side of the 
curves is excluded. 
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4. The transition probability e(z/^ z/g) as a function of the detection 
threshold energy Ett of muons for sin^ 29 = 0.5 and A/ = 0, which is 
solely due to the standard MSW effect. The solid and dashed curves 
have parameters Am^ = 1 x 10~^eV^ and 0.5 x 10~^eV^, respectively. 

5. The transition probability e(i/^ — > fg) as a function of the detection 
threshold energy E'th of muons for sin^ 29 — 0.5 and A/ = —0.5 x 10"^^. 
The solid, dashed and dotted curves have parameters Am^ = 1 x 
lO-^eV^ 0.5 X 10-2eV^ and 0, respectively. 
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